Renormalization of the Deuteron with One Pion Exchange 
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We analyze the deuteron bound state through the One Pion Exchange Potential. We pay at- 
tention to the short distance peculiar singularity structure of the bound state wave functions in 
coordinate space and the elimination of short distance ambiguities by selecting the regular solution 
at the origin. We determine the so far elusive amplitude of the converging exponential solutions at 
the origin. All bound state deuteron properties can then be uniquely deduced from the deuteron 
binding energy, the pion-nucleon coupling constant and pion mass. This generates correlations 
among deuteron properties. Scattering phase shifts and low energy parameters in the 3 Si — 3 Di 
channel are constructed by requiring orthogonality of the positive energy states to the deuteron 
bound state, yielding an energy independent combination of boundary conditions. We also analyze 
from the viewpoint of short distance boundary conditions the weak binding regime on the light of 
long distance perturbation theory and discuss the approach to the chiral limit. 
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I. INTRODUCTION 



Pion dynamics plays a dominant role in the low en- 
ergy structure of the Nucleon-Nucleon interaction, and 
in particular in the description of light nuclei like the 
deuteron yj. The long distance part of the interaction 
is given by one, two and higher pion exchanges and the 
fact that the deuteron is a weakly bound state suggests 
that many of its properties can indeed be explained in 
terms of these dynamical degrees of freedom in a model 
independent way and regardless on the less known short 
distance interaction. Glendenning and Kramer [2] in 
the early sixties recognized clear correlations between 
several deuteron observables generated by truncating 
the One Pion Exchange (OPE) potential at a distance 
R = 0.4915fm and assuming a hard core inside. Tight 
constraints on deuteron observables were established by 
Klarsfeld, Martorell and Sprung 0, by integrating the 
deuteron wave function from infinity down to a cut-off 
radius using the OPE potential and rigorous inequali- 
ties. An accurate determination of the D/S asymptotic 
ratio was made by Ericson and Rosa-Clot 0, 0| based 
on assuming the OPE correlation between the S and D 
wave functions and taking realistic potential models to 
describe the S wave function, (for a review on these 
developments see ej* Ref. [7j). Friar et al. use a mul- 
tipole form factor [8j whereas Ballot et al. used separate 
monopole forms factor for the central and tensor part 
of the OPE potential mimicking the finite size of nucle- 
ons 0, 0| . Along a similar line of investigation Sprung 
et al. used a square well potential yj for the central 
component and a vanishing potential for the tensor com- 
ponent. 
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Within the effective field theory (E FT) approach to 
nuclear physics proposed by Weinberg [12J the situation 
was revisited from a somewhat different perspective since 
the OPE potential appears as the lowest order of a per- 
turbative hierarchy based on chiral symmetry [l3| (for a 
review see e.g. Ref. Q5), and short distance ambiguities 
could be eliminated by the renormalization program if 
the auxiliary regulator is removed from the theory at the 
end of the calculation. This cut-off independence should 
occur at any level of approximation, no matter how many 
pions are exchanged. At long distances, renormalization 
group methods suggest that one is close to an infrared 
fixed point |15| . The renormalization procedure can be 
explicitly and analytically carried out within perturba- 
tion theory [lq|. However, these nice features become a 
non trivial numerical problem beyond perturbation the- 
ory motivating the use of truncation cut-off schemes. The 
work of Ref. [TtJ uses a gaussian cut-off in coordinate 
space to regulate the contact delta interaction, Ref. y~3 
proposes the use of a subtraction method in momentum 
space regulating the central part, Ref. yjj uses a sharp 
momentum cut-off and in Ref. _2Qj it was pro pos ed to 
use a finite short distance cut-off, whereas Ref. 21] puts 
exponentially suppressed regulators in momentum space. 
It should be mentioned that in all cases the correspond- 
ing coordinate/momentum space cut-off parameter a/ A 
is uncomfortably large/small from the viewpoint of renor- 
malization theory. Typically, one has a ~ 1.4fm (see e.g. 
Ref. |3) and A = 600MeV yj| respectively. So, it is not 
obvious that according to the basic principles of EFT the 
short distance ambiguities are, as one might expect, in- 
deed under control. Moreover, the existence of a well 
behaved finite renormalized limit is never guaranteed a 
priori and one relies mainly on numerics. Actually, the 
fact that the results on deuteron observables look rather 
similar, regardless on the particular way how the poten- 
tial is modeled at short distances, proves that the long 
distance pion dynamics dominates the physics confirming 
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the findings of Glendenning and Kramer 2] 

more than 40 

years ago, but does not resolve the mathematical problem 
whether the OPE potential can make unambiguous pre- 
dictions regardless of any short distance physical scale. 

The OPE potential is local in coordinate space where 
the problem is naturally formulated by the standard 
Schrodinger framework. Moreover, it is singular at the 
origin and giving boundary conditions at that point is 
not a well defined procedure for uniquely determining 
the solution |22j ] (for a comprehensive review in the one 
channel case see e.g. Ref. [23] )■ There is the added 
difficulty that we have two coupled second order differ- 
ential equations. In the deuteron channel one has four 
independent solutions, which according to their singular- 
ity structure correspond to cither two regular and two 
irregular solutions at infinity or three regular and one 
irregular solution at the origin. The normalizability con- 
dition of the deuteron wave functions eliminates all con- 
stants for a given deuteron binding energy, which in- 
stead of being predicted has to be treated as an inde- 
pendent parameter. One of the advantages of the coor- 
dinate space treatment of renormalization is that it can 
directly be extended to other singular cases such as the 
Two Pion Exchange (TPE) potential which is also 
finite everywhere except the origin. In contrast, momen- 
tum space treatments require an extra regularization of 
the potential besides the standard cut-off regularization 
of the Lippmann-Schwinger equation. 

The authors of Ref. 01 found a discrete sequence of 
equivalent short distance cut-off radii having almost the 
same deuteron properties. In their analysis of the prob- 
lem one regular solution at the origin with a converging 
exponential behavior, exp(— 4(2i?/r)2 ) with R ~ lfm, 
was discarded on numerical grounds. The same result 
was also implicitly used in Ref. |2^| and larg e N c argu- 
ments in favor of it were raised in Ref. 26] . This ex- 
tra condition would actually predict the deuteron bind- 
ing energy from the OPE potential. As we will show in 
this paper, the converging exponential is non-vanishing 
although rather elusive because its contribution to the 
deuteron wave function only becomes sizeable at rela- 
tively large distances and accurate numerical work must 
be done to pin down its value with a certain degree of 
confidence. 

In the present work we show that there is no need to 
truncate the OPE potential on a physical scale to produce 
unique and cut-off independent predictions for deuteron 
properties and scattering observables in terms of the OPE 
potential parameters and the deuteron binding energy. 
These might then legitimately be called OPE model in- 
dependent predictions and paves the way for a systematic 
investigation on the case where more pions are exchanged 
and other effects are taken into account |24|. After pre- 
senting the basic notation in Sect. Ill Al we discuss in 
Sect. Ill Bl the regular solutions at the origin and estab- 
lish that the limit when the regulator is removed is finite. 
For numerical purposes it is useful to define some short 
distance regulator as an auxiliary tool. In Sect. Ill CI we 



use six different regulators based on boundary conditions 
of the wave function and check for stability to high preci- 
sion for all regulators. This procedure generates correla- 
tions among deuteron observables if the deuteron binding 
energy is varied as a free and independent parameter as 
we do in Sect. ITTTl A particularly interesting situation is 
provided by the weak binding limit, which can be taken 
with fixed OPE potential parameters. In such a case the 
long distance behavior should dominate and one might 
expect perturbative methods to apply and be compared 
to the exact OPE calculation. The details of the pertur- 
bative calculation are postponed to Appendix^Jwhere we 
present a coordinate space version of the method, in con- 
sonance with the exact treatment. A detailed comparison 
shows that the perturbative argument is too naive and 
would only hold in the weak coupling regime as well, due 
to the appearance of non-analytical contributions in the 
ttNN coupling constant. Mathematically, we show that 
it is not possible to go beyond first order since the co- 
efficients of the expansion diverge. Numerically, the dis- 
agreement at first order is typically on the 30% level for 
physical values of the OPE parameters at zero binding. 
After Ref. [2]J the chiral limit in nuclear physics has at- 
tracted considerable attention in recent works [28Ll29H30| 
and also the limit of heavy pions in connection with lat- 
tice QCD calculations, where the pion mass is still far 
from its physical value. We study in Scct. lIVI the correla- 
tions among those observables if the pion mass is varied 
away from its physical value by studying a suitable exten- 
sion of the Feynman-Helmann theorem. Another remark- 
able property of the OPE potential which we deal with 
in Sect. [V]is that low energy parameters as well as the 
scattering phase-shifts can be uniquely determined from 
the OPE potential parameters and the deuteron binding, 
due to orthogonality constraints of the bound state and 
scattering states. In Sect. IVII the determination of the 
non-vanishing coefficient of the converging exponential 
at the origin is carried out by a short distance expansion 
to eighth order of the OPE deuteron wave functions. In 
Sect. IVIll we come to the conclusions. 

One of the surprising results in the OPE description of 
the deuteron has to do with the small asymptotic ratio 
between the D and S waves, w(oo)/u(oo) = r) = 0.0256 
coming from a large ratio at short distances of order 
unity, w(0)/u(0) = l/y/2 = 0.707. Although this fea- 
ture is specific to the OPE potential it is somewhat a 
bit outside the main topic of this work. So we relegate 
this issue to Appendix ICl where we show how this can be 
easily understood if a local rotation of the deuteron wave 
functions diagonalizing the coupled channel potential is 
carried out. Obviously, such a transformation cannot si- 
multaneously diagonalize the kinetic terms, but the resid- 
ual mixing is related to the derivative of a local mixing 
angle which numerically turns out to be a slowly varying 
function. Using this as a starting approximation we can 
determine in a perturbative fashion the asymptotic D/S 
ratio yielding the exact OPE value with a 1% accuracy. 
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II. BOUND STATE EQUATIONS AND THEIR 
SOLUTIONS 



A. The OPE deuteron Equations 

The Deuteron coupled channel 3 Si — 3 Di set of equa- 
tions read 

- u"(r) + U s (r)u(r) + U s d(r)w(r) = -^ 2 u{r) , 

(1) 



-w"(r) + Usd{r)u{r) + 



U d (r) 



w(r) = —^ 2 w{r) , 



(2) 



together with the asymptotic conditions at infinity 
u(r) - A s e-v, 
w(r) -> A D e-< r ( 1 + — + 



(3) 



where 7 = v 'MB is the deuteron wave number, As is 
the normalization factor and the asymptotic D/S ratio 
parameter is defined by 77 = Ap/Ag- The 3 Si — 3 Di 
coupled channel potential is given by 



U sd = 2V2U 7 



U d = Uc-2U T , (4) 



where the OPE reduced potential (U — 2fiV) is given for 
r > by 



U c 
U T 



mMg 2 



16tt/2 r 
m 2 Mg 2 A e 
~^T 2 



3 3 
1 + — + 7 — ^ 



(5) 
(6) 



where m is the pion mass, M = 2/i np = 2M„M P / (M n + 
M p ) twice the reduced proton-neutron mass, gA the axial 
nucleon coupling constant and the pion weak decay 
constant. Note that we assume this potential to be valid 
for any strictly positive distance, r ^ 0, so the limit 
r — > + will be carefully taken, without subtracting any 
contribution in the potential. 

It is convenient to define the length scale 



R 



327r/-? 



(7) 



which value is around 1 fm. For numerical calculations we 
take f n = 92.4MeV, M = 938.918MeV, g A = 1.29 and 
hence R = 1.07764fm. Using the Goldberger-Treiman 
relation, g^NN = 9AM/fn, the corresponding pion nu- 
cleon coupling constant is g^NN = 13.1083 according to 
a phase shift analysis of NN scattering jsj. Neverthe- 
less, after the latest determinations from the GMO sum 
rule we will also take the value g„NN — 13.3158. As 
we will see this variation at the 5% level dominates the 
uncertainties in the OPE calculations. 



B. The short distance regular solutions 



We look for normalized functions of the Eqs. J2J, 



1 = 



(u(r) 2 + w{r) 2 )dr , 



(8) 



from which As can be determined. The normalization at 
infinity is guaranteed due to the asymptotic conditions, 
Eq. |J3}. However, the coupled channel potential becomes 
singular at short distances, since Ut — > —2R/r 3 . Keep- 
ing only this term in Eqs. J2J) one can decouple the equa- 
tions through the unitary transformation |ll| 



u A (r) 
ur{t) 



12 1 

3«(r) + ^f"( r J ' 



1 



V3 



u(r) 



w(r) , 



(9) 



yielding an attractive singular potential Ua — > — 4i?/r 3 
for ua and Ur — > 8R/r 3 for ur. Any solution obtained 
by integrating from infinity with the Eq. (j3J) down to the 
origin has the asymptotic short distance behavior , 



(10) 



u A (r) 



r \ 3/4 r 
R 



C 1A e-^- +C 2A e 



The constants C\r, C2R, C\a and C2A depend on both 
7 and 77 and the OPE potential parameters, g-nNN and 
to. Note that the leading short distance r dependence 
does not involve the pion mass and the deuteron wave 
number. Higher order corrections to these solutions can 
be computed systematically to high orders and are pre- 
sented below in Sect. I VII 

The regular solution at infinity contains the normal- 
ization constant As, which is customarily set to one for 
computational purposes, the deuteron wave number 7 
and the asymptotic D/S ratio parameter 77. The normal- 
izability of the wave function at the origin requires 



(11) 



which is a relation between rj and 7. The other remain- 
ing constants are then completely fixed. This means that 
for the OPE potential, the deuteron binding energy can 
be used as an independent parameter. Thus, one has 
three independent variables, 7, the coupling constant 
with length scale dimension R (or equivalently g n NN) 
and the pion mass m. Obviously, this suggests integrat- 
ing in from infinity and determining 77 from the regularity 
condition at the origin 1)110. 



The solutions fo r u a and uji are written in terms of spherical 
Bessel functions We keep the leading short distance behav- 
ior only. 
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To analyze whether some additional condition arises 
let us check the selfadjointness of the coupled channel 
Hamiltonian. The flux at a point r is given by 



iJ(r) = u* (r) u(r) — u* (r)u (r) 

+ w* (r)'w(r) — w* (r)w' (r) , 



(12) 



so that current probability conservation at the origin im- 
plies 



\C 



IA 



\C 2 a\ 2 = 2V2i (C* 1R C 2R - C* 2R C 1R ) . (13) 



Thus, if we set C\ R = there is no condition on C 2R and 
one has C\a — Ca^ w and C 2 a — CAe~~ lip with Ca and 
ip real. So, we have three constants, C 2R (j), Ca{i) and 
95(7), characterizing the normalizable solutions at short 
distances for a given value of the deuteron wave number 

7, 



u R {r) 
u A {r) 



C, 



,(7) (£ 



x 3/4 
) ' 



-4V2V- 



3/4 



sin 



(14) 



Actually, if we have any other state, say a scattering 
state with positive energy, unitarity (i.e. orthogonality) 
requires that the constant ip(k) coincides with the bound 
state phase ^(7). We will come back to this issue later 
when discussing low energy parameters and scattering 
solutions in Sect. It is natural to expect that some 
combination of short distance constants is independent 
on the OPE potential parameters as they encode short 
distance physics. In Sect. lIVl we establish, by demanding 
the standard Feynman-Hellmann theorem, that specifi- 
cally the short distance phase cp does not depend on the 
OPE potential parameters. In Sect. lVII we determine the 
values of the three constants characterizing the three reg- 
ular solutions by a detailed short distance analysis of the 
OPE deuteron wave functions. 

Note that any additional condition would actually pre- 
dict both 7 and 77 from m and R. This contradicts the 
claim of Ref . fill that C 2R = 0, a conclusion implicitly 
used in Ref l25l and supported by the large iV c argu- 
ment of Ref. |2g. On the other hand, if one takes the 
experimental values of r\ and 7 as done in Ref. [2(j one 
obtains both non vanishing C\ R and C 2R i.e., the irreg- 
ular non-normalizable solution, unless a short distance 
cut-off, R > 0.8fm, is introduced as a physical scale and 
not as an auxiliary removable regulator. 



run into the diverging exponentials due to round-off er- 
rors for r ~ 0.05fm and dominate over the converging 
exponential. The reason has to do with the fact that the 
natural scale where both exponentials are comparable is 
rather large r = 4^/2R ~ 6fm, but in that region the 
lowest order short distance approximation does not hold. 

Instead, we will also try putting several short distance 
boundary conditions corresponding to the choice of reg- 
ular solutions at the origin, 

u(a) = (BC1), 
u'(a) = (BC2), 
w(a) = (BC3), 
w'{a) = (BC4), 
u(a )-V2w(a) = (BC5), 
u'(a) — V2w'(a) = (BC6), 

(15) 

The advantage of using this kind of short distance cut- 
offs based on a boundary condition is that there is only a 
single scale in the problem as one naturally expects, and 
that one never needs to declare what is the wave function 
below the boundary radius. Putting a square well poten- 
tial as a counter-term [25| with depth Uq appears natural 
from standard perturbative experience but needs speci- 
fication of a further length scale, 1/^/Uq, and moreover, 
generates multi- valuation ambiguities |3J, l35| . 

It is convenient to use the superposition principle of 
boundary conditions to write 



u(r) = us{r) + -quD{r) 
w(r) = u>s{r) + rjWD{r) , 



(16) 



where (us, ws) and (ud, wd) correspond to the boundary 
conditions at infinity, Eq. (j3Jl with Ag — 1 and A^, = 
and with As = and Ap = 1 respectively. Thus, at the 
boundary we can impose any of the conditions by just 
eliminating 7] 2 . The resulting 77 value obtained by all 
these boundary conditions is presented in Fig. ^ Actu- 
ally, we see that the boundary condition u'(a) — \/2w' (a) 
is about the smoothest condition we can think of, since 
the u R combination goes to zero at small distances, its 
derivative, u' R also goes to zero, although a bit less faster 
since u' R /u R ~ 1/r 3 / 2 . We see that all determinations of 
r] based on any of the proposed cut-offs yield the same 
value with great accuracy at cut-off radii below 0.2fm. 



C. Regularization with boundary conditions 

Ideally, one would integrate in the large asymptotic so- 
lutions, Eq. ©, and match the short distance behavior 
of Eq. I|llfl imposing the regularity condition l|ll|) . In 
practice, however, the converging exponential at the ori- 
gin is rather elusive since integrated-in solutions quickly 



2 Numerically we find at the cut-off boundary r = 0.2fm 

u(0.2) = 1139.23 -43263.2 r) 

w(0.2) = -1807.33 + 68632.5 r\ 

m'(0.2) = -35529.8 + 1.34913 ■ 10 6 »? 

u/(0.2) = 55194.3 - 2.09606 ■ 10 6 ?7 (17) 

These large numbers appear because the of the dominance of the 
diverging exponential at short distances. 
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An important observable is the deuteron inverse radius 



0.15 0.2 0.25 0.3 0.35 0.4 
Cutoff Radius a (fm) 



0.45 0.5 



FIG. 1: The dependence of the asymptotic D/S ratio r\ on the 
cut-off radius a for several boundary conditions. We use m = 
138.03MeV and R = 1.07764fm (corresponding to g nNN = 
13.1083). 



This is somewhat fortunate since arithmetic precision is 
outraged typically for r < 0.06fm. Obviously, any short 
distance cut-off generates finite cut-off effects in the wave 
functions for distances close to the cut-off radius. In 
Sect. IVIl we analyze this problem by matching solutions 
of the form of Eq. I|14|) to the integrated in numerical 
solutions and find that for many practical purposes these 
finite cut-off effects are negligible. Thus, we will base 
most of our results on the "smoothest" condition BC6 of 
Eq. O. 

The resulting deuteron wave functions u and w ob- 
tained by integrating in from infinity to the origin the 
OPE potential are plotted in Fig. [3 where the irregular 
solutions are obtained with the experimental D/S ratio 
r]d = 0.0256 and the regular ones with the OPE D/S 
ratio 7/ope = 0.026333. For comparison we also plot the 
Nijmll deuteron wave functions. We emphasize that the 
value of rj is a direct consequence of taking the OPE down 
to the origin seriously. 



D. Deuteron observables 

Once the solutions are known we can determine several 
observables of interest. The matter radius reads, 



1 1 r°° 

r ™ = \^ = \\ r 2 (u(rf+w(rf)dr (18) 
•J o 

while potential contribution to the quadrupolc moment 
(without meson exchange currents) 



i r°° 

?d= 20/ r2w ( r )( 2 ^ u ( r ) - Mr))dr (19) 



dr 



u(r) 2 + w(r) 2 



(20) 



which appears in low energy pion-deuteron scattering. 
Finally, the -D-state probability is given by 



Pd = 



w(r) 2 dr 



(21) 



Both Pd and (r _1 ) are sensitive to the intermediate dis- 
tance region around 2fm whereas Qd and r m get their 
contribution from larger distances ~ 4fm. 

The results for the asymptotic S-wave normalization 
As, the matter radius r m , the quadrupole moment, Qd, 
and the D-state probability, Pd are presented in Table [I] 
The errors in the numerical calculation have been as- 
sessed by varying the short distance cut-off in the range 
a = 0.1 — 0.2fm (in momentum space that would naively 
correspond to take A=l/a = 2 — 4GeV) . As we see, 
the cut-off uncertainty is smaller than the one induced 
by variations at the 2% level in the g^NN coupling con- 
stant in the range between the lowest value (~ 13.1) ob- 
tained by a fit to NN phase-shifts [3l| and the highest 
recent value (~ 13.3) determined from the GMO sum 
rule |32] ]. Equivalently, this uncertainty corresponds to 
take R — 1.0776fm and R = 1.1108fm respectively. Our 
results are generally speaking in agreement with previous 
determinations where different sorts of cut-off methods 
have also been implemented. 



E. Discussion 

At this point it may prove useful to ponder on the pre- 
vious results from a wider perspective. Let us remind 
that the basic assumption of an EFT is that the study 
of long wave length phenomena such as low energy scat- 
tering or weakly bound systems do not require a detailed 
knowledge of short distance physics. This general and 
widely accepted principle requires some qualification be- 
cause attractive and repulsive singular potentials behave 
quite differently in this respect. Singular attractive po- 
tentials, ~ l/r n , generate wave functions vanishing as a 
power law, r n / 4 sin(r~™/ 2+1 +ip), and which need a mixed 
boundary condition to specify the short distance phase 
(p. Thus, short distance details become less important, 
regardless on the value of (p. On the contrary, for sin- 
gular repulsive potentials the wave functions behave as 
r n/4 e ±r 1 + anc j on iy f or t ne regular solution short dis- 
tance details become irrelevant. In the OPE potential, 
it is precisely the repulsive short distance OPE compo- 
nent which requires a fine tuning of the solutions and 
eliminates one a priori independent parameter like, e.g., 
the asymptotic D/S ratio r\. As we see, if 77 is treated 
as an independent variable the short distance behavior 
of the deuteron wave functions precludes the definition 
of a normalizable state due to the onset of the irregular 
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-0.2 ' Deuteron Wave Functions 

I , i i , 



1 2 3 4 5 

r(fm) 

FIG. 2: The deuteron wave functions u and w, obtained by 
integrating in from infinity to the origin the OPE potential, 
compared to those obtained with the Nijm II potential [33|. 
The irregular solutions are obtained with the experimental 
D/S ratio, r\d = 0.0256, and the regular ones with the OPE 
D/S ratio, ijope = 0.026333. We use m = 138.03MeV and 
R = 1.07764fm (corresponding to q^nn = 13.1083). 

solution. This short distance insensitivity at low energies 
could only be implemented by keeping the experimental 
77 value and ignoring OPE physics below some scale. The 
lower limit established in Ref. [2(J to obtain a normaliz- 
able state was a ~ 1.3fm for the OPE potential. This 
obviously requires some extension of the wave function 
below that scale and the pretended model independence 
becomes a bit obscured. Our point is that the short dis- 
tance insensitivity materializes automatically for the reg- 
ular OPE deuteron wave functions since they vanish at 
the origin. 



III. OPE CORRELATIONS IN DEUTERON 
OBSERVABLES 

As we have said, in the OPE potential we can use the 
deuteron wave number as an input of the calculation on 
the same footing as g^NN arL d the pion mass m. Then, 
other observables are predicted. We will study now the 
dependence of these observables on 7, m and R. 

A. Dependence on the Binding Energy 

In Fig. |31 we show the dependence of the D/S ratio 
as a function of the deuteron wave number 7 keeping 
to and R fixed. In the weak binding limit 7 <C m^, 
long distances dominate and the finiteness of the wave 
function at a point r 3> 1/m requires 77 ~ j 2 . The ra- 
dius of convergence of such an expansion for the observ- 
ables is I7I < to/2, since the integrals involve the factor 



e -(27+m)r ^ large distances, diverging for 7 < —to/2. 
The experimental number is not far from 7 = m/3, which 
is within the domain of analyticity but somewhat close 
to the convergence radius. So, one expects a slow con- 
vergence. As we see in the weak binding limit we have a 
quadratic behavior ?7ope ~ 7 2 whereas for stronger bind- 
ing a linear behavior sets in. It is remarkable that the 
experimental values in the intermediate regime. On the 
other hand, in the strong binding case 7 S> m^, short 
distances dominate and we must have r\ ~ 1/V2. Nu- 
merically we find for the deuteron observables, 

r/ OPE = 0.96387 2 - 3.46864 7 3 + C( 7 4 ) (22) 
4OPE 

-f=- = 1 + 1.24557 - 0.4705 7 2 + 0( 7 3 ) (23) 
V 2 7 

^7 r° PE = 1 + 1-24557 ~ 0.47057 2 + C( 7 3 ) (24) 

QO pE = 0.6815- 3.54377 + C(7 2 ) (25) 

Note that we have the weak binding correlation 

r m = ^§j, + 0(7 3 ) (26) 

which is compatible at the 2a confidence level with data; 
for the experimental value As — 0.8845(8) the value 
r m = 1.984(2) to be compared with the experimental 
number r m = 1.971(6). In the weak binding limit one 
also has the correlation 

^^ = 1 + 0(7) (27) 
fid 

a dependence that one would expect on general grounds 
by just taking the asymptotic formulas and neglecting the 
w(r) 2 term in the expression for the quadrupole moment. 
Experimentally this relation is fulfilled with a 15% accu- 
racy. The non-perturbative OPE value is actually closer 
to potential models. 



B. Comparison with perturbation theory 

It is instructive to solve the coupled deuteron equa- 
tions, Eq. Q in standard perturbation theory for the 
fixed energy bound state. One of the reasons is to check 
the correctness of our non-perturbative calculations in 
the weak binding regime. Another motivation is to estab- 
lish contact with the perturbative calculations of Ref. 0] 
where dimensional regularization in the power divergence 
subtraction (PDS) scheme was implemented. Finally, 
there is the question of quantitatively assessing the valid- 
ity of such an approximation. We relegate the calculation 
to Appendix IaI At first order in perturbation theory one 
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TABLE I: Deuteron properties for the OPE compared to the short range approximation and first order perturbation theory. We 
use the non-relativistic relation 7 = ^/2u np B with B = 2.224575(9) and take m = 138.03MeV and R = (3/8M)gl NN /(4vr) = 
1.07764fm corresponding to g^NN ~ 13.1083 |3 111 (except the row OPE* where the value g^NN = 13.316 [32| has been taken). 
The error is estimated by changing the short distance cut-off in the range a = 0.1 — 0.2fm. 





7 (fm" 1 ) 


V 


A s (fm 1/2 ) 


r m (fm) 


Qd(fm 2 ) 


Pd 


(r J ) 


a (fm) 


a 02 (fm 3 ) 


a 2 (fm 5 ) 


r (fm) 


Short 


Input 





0.6806 


1.5265 





0% 


00 


4.3177 











OPE(pert) 


Input 


0.051 


0.7373 


1.6429 


0.4555 


0% 


00 


4.6089 


2.5365 





0.4831 


OPE 


Input 


0.02633 


0.8681(1) 


1.9351(5) 


0.2762(1) 


7.88(1)% 


0.476(3) 


5.335(1) 


1.673(1) 


6.169(1) 


1.638(1) 


OPE* 


Input 


0.02687 


0.8718(2) 


1.9429(6) 


0.2826(2) 


7.42(1)% 


0.471(3) 


5.353(1) 


1.715(1) 


6.4001(1) 


1.663(1) 


Nijmll 


Input 


0.02521 


0.8845(8) 


1.9675 


0.2707 


5.635% 


0.4502 


5.418 


1.647 


6.505 


1.753 


Reid93 


Input 


0.02514 


0.8845(8) 


1.9686 


0.2703 


5.699% 


0.4515 


5.422 


1.645 


6.453 


1.755 


Exp. a 


0.231605 


0.0256(4) 


0.8846(9) 


1.9754(9) 


0.2859(3) 


5.67(4) 




5.419(7) 






1.753(8) 



a (Non relativistic) . See e.g. Ref. LSTjl and references therein. 



gets in the weak binding limit 



?7pcrt 
^S,pcrt 

r - m , P ertV / 87 

Qpcrt 



1.54977 2 - 4.154797 3 + C(7 4 , R 2 ) 
1 - 0.7184 7 - 2.7394 7 2 + C( 7 3 , R 2 ) 

1 + 0.71843 7 - 2.7394 7 2 + 0{~/ 3 , R 2 ) 
1.09587 - 5.87576 7 + 0{"f 2 ,R 2 ) 

(28) 



The nominally 0{R 2 ) second order contributions are in 
fact divergent because the leading order correction to the 
D-wave component w(r) diverges at the origin (see Ap- 
pendix^). In general terms we find that the exact OPE 
results are estimated within 30% by the first order per- 
turbative calculations of Appendix. lAl 



DEPENDENCE ON THE PION MASS AND 
CHIRAL LIMIT 



IV. 



Recent works |27J, |2J, |2<J, |30j predict the change of the 
deuteron binding energy as a function of the pion mass 
by taking the experimental binding energy at the phys- 
ical value of the pion mass and making the additional 
assumption short distance physics to be independent on 
the pion mass. While it is true that the leading short 
distance r dependence of the deuteron wave functions 
are independent on the pion mass, the constants Ca, Cr 
and ip do in principle depend on the three independent 
parameters to, g^NN and 7. As we have noted 7 can- 
not be predicted for the OPE potential. So the approach 
pursued in Refs. [2?i l2ci l2fjl l30| is equivalent to inte- 
grate in with the physical pion mass and then integrate 
out fixing some combination of short distance constants 
with the unphysical pion mass and searching for the ap- 
propriate regular solution at infinity. If one makes the 
pion lighter long distance effects should dominate, and 
one could just use the OPE potential to estimate the 
chiral limit as a first approximation. It is thus interest- 
ing to analyze the pion mass dependence both explicitly 



(i.e. varying to in the OPE potential) and implicitly (i.e. 
taking into account the dependence of the OPE coupling 
R on the pions mass). We will determine the pertinent 
combination of short distance constants by demanding 
the Feynman-Hellmann theorem in the OPE potential. 



A. Explicit pion mass dependence 

To proceed, let us assume that to an infinitesimal 
change in m — > to + Am there corresponds a change 
both in the deuteron wave number 7 — > 7+ A7 and in the 
coupled channel potential matrix U(r) — > U(r) + AU(r). 
We can write a Lagrange identity by varying the equation 
and its adjoint, yielding for a normalized state 



dm 



+ 



dU , 
dm 
, du 
dm dm 



du' 



, dw 
dm 



dw' 
dm 



(29) 



This is an extended Feynman-Hellmann theorem where 
the second term in the l.h.s. corresponds to the short 
distance contribution (the term at infinity vanishes for 
a bound state). One of the advantages of the Feynman- 
Hellmann theorem is that one could in principle establish 
comparison theorems, provided the change in the cou- 
pled channel potential matrix, AU, is a definite quadratic 
form. Note also that the derivative with respect to to an- 
nihilates the centrifugal term, 6/r 2 , and one can diago- 
nalize the coupled channel potential by the unitary trans- 
formation Eq. (O so that the result behaves additively 
in the attractive and repulsive eigenchannels. Using the 
leading short distance behavior, Eq. I|14[l . we therefore 
get 



dj^ 
dm 



dr 







u A {r) 



.dU A 
dm 



un{r) 



,dU E 



dm 



C 2 dl P 
A dm 



(30) 
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FIG. 3: The dependence of the asymptotic D/S ratio r\ (upper left panel), the S-wave normalization As (in fen -1 / 2 , upper 
right panel) , the matter radius r m (in fm, lower left panel) and the quadrupole moment Qd (in fm 2 , lower right panel) on the 
deuteron wave number 7 (in fm - ) for the short range theory, the first order perturbative result and the exact OPE result. 
The points represent the experimental values. We use m = 138.03MeV and R — 1.07764fm. 



As we see, assuming as suggested in Ref. [23 that the 
short distance physics does not depend on the pion 
mass corresponds to demanding the standard Feynman- 
Hellmann theorem where only the OPE potential change 
contributes. For Ca ^ one obtains the condition 



d dip g?7 ^ dip 

dm ' dj dm dm 



0, 



(31) 



whence a functional relation between the pion mass and 
the deuteron binding energy follows. However, note that 
even in this case the sign of the result is indefinite since 



dm 



dr 



u A (r) 



2 dU A 
dm 



ur{t) 



,dU h 



di 



(32) 



So, we have to determine the sign numerically. 

The relation (|31|l has an equivalent formulation in the 
boundary condition regularization. For instance, if we 
assume the same condition BC6 of Eq. (|43|) for all values 
of the pion mass we get 



_d_ 

dm 



w'{a) 



u(a)^/2 + uu(a) 



= 0. 



(33) 



where a is taken to be independent of m. In practice, one 
computes the ratio within the bracket for the physical 
pion mass and searches for 7 such that the ratio for the 
unphysical pion mass yields the same numerical value. If 
we take the chiral limit we get Bd(0, g^NN) = 4.3539MeV 
This value is very close to the one found in Ref. [2{| 
Bd(0,g7rNN) = 4.2MeV 3 . Deuteron observables in the 
explicit m = limit are listed in Table ITT1 



3 These authors look for poles of the S— matrix so constructed 
as to reproduce the physical scattering length oq = 5.42fm and 
effective range ro = 1.75fm at the physical value of the pion mass 
in the 3 Si eigen channel. By doing so the explicit dependence on 
9ttNN becomes rather weak. Actually, they take g-^NN = 12.73 
and we would get instead -8^(0, g^NN) = 0.98MeV instead. This 
apparent contradiction is resolved by noting that, as we will see 
below, for g n NN = 13.1083 in the OPE we get an scattering 
length 00 = 5.335fm an effective range of ro = 1.63 quite close 
to the experimental values. 
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B. Implicit pion mass dependence 



To take into account the implicit pion mass depen- 
dence we have to take into account the dependence of 
R = 3g A AI/32irf% on the pion mass. In the chiral limit 
one gets a larger OPE coupling |29j ■ The value is uncer- 
tain and as an educated guess we take Rq = 1.06(2)i?. 
Using the same formulation as in the m dependence, the 
change in the deuteron binding with respect to the g n NN 
coupling constant or equivalently the scale dimension R 
we get (assuming as before the short distance angle (p to 
be independent on m), 




20 40 60 80 100 120 140 
Pion Mass ( MeV) 



R ^R = I dr[u A (r) 2 U A + u R (r) 2 U B ] 



(34) 



FIG. 4: Explicit pion mass dependence of the deuteron 
binding energy, keeping g n NN = 13.1083 fixed. We have 
B d {m) = 2.22457MeV, B d (0) = 4.3539MeV. 



Again, the result is indefinite since Ua < and Ur > 
and it is not obvious, unlike naive expectations, that a 
stronger coupling provides stronger binding. The sign 
depends actually on the details of the wave functions and 
the particular values of the parameters. Numerically one 
finds drf/dR > 0, a trend that can be understood if the 
repulsive term in Eq. l|34[l is neglected or on the basis of 
the inequality \ua\ > \ur\ which is numerically fulfilled. 
In any case one has the differential inequality dj/dR < 
j/R. 

Numerically, we get 70 = 0.61(10)fm~ 1 and hence 



B Q d = 15(5)MeV 



(35) 



a value compatible with the analysis of Ref. [23 = 
9.6 ±3 (perhaps with larger errors [30l | 4 ) but in disagree- 
ment with Refs. |25ll28l] where the deuteron becomes un- 
bound for to < 90MeV. In any case we confirm the trend 
of having a stronger binding of the deuteron in the chiral 
limit. The corresponding observables can be looked up 
in Table. M 



SCATTERING PROPERTIES IN THE 
3 Si - 3 £>i CHANNEL 



A. Orthogonality constraints and Phase Shifts 

For the a and (3 positive energy scattering states we 
choose the asymptotic normalization 



Uk, a (r) 

Wk,a{r) 

Wk,p(r) 



-r-r(jo(kr) cos<5i -y (kr) sin Si 
smd-y V 

-\h(kr) - mikr) sin (5] 
1 v 



sin 
sin e 
sin S 



(36) 



1 



sin Si 
tane 
sin<5i 



jo (At) cos<5 2 - ?/o (At) sin £ 



j 2 (At) cos 8 2 - Vi (At) sin 5 ; 



(37) 



4 If we take Ro 



1.1R we get B° d = 33MeV. 



where ji (x) — xji (x) and yi (x) = xyi (x) are the re- 
duced spherical Bessel functions and Si and 82 are the 
eigen-phases in the 3 Si and 3 Di channels, and e is the 
mixing angle E%. Again, the general solution at short 
distances is given by the general Eq. I|14|l , where the con- 
stants Ca-, Cr and ip are now different since we have 
a zero energy state and depend whether we have an 
a or (3 state, so we have the short distance constants, 
CA,a(k), C R)0l (k), ip a {k) and CA,p(k),C Rt p(k),<pp(k) re- 
spectively. This implies certain correlations between 8±, 
82 and e. 

For a regular self-adjoint potential the orthogonality of 
bound and scattering states comes out automatically. We 
look now for the consequences of demanding this prop- 
erty in the singular OPE potential. Using the standard 
manipulations to prove orthogonality between states of 
different energy we get the following relation between a 
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TABLE II: Deuteron properties for the OPE and their dependence on the pion mass. We use the non-relativistic relation 
7 = ^2yL nv B = 0.231605 far 1 with B = 2.224575(9) for m = 138.03 and g nNN = 13.1083. m = (explicit) means taking 
m = but using R — 1.07764fm (or equivalently g-nNN = 13.1083). m = (implicit) means taking m = and using 
Ro = 1.06(2) x R. 





7 (fm- 1 ) 


V 




r d (fai) 


Qd(fm 2 ) 


Pd 


m = 138.03 MeV 
m = (explicit) 
m — (implicit) 


Input 
0.3240(1) 
0.61(10) 


0.02633 
0.09452 
0.15(2) 


0.8681(1) 
0.8444(1) 
0.48(7) 


1.9351(5) 
1.550(1) 
0.98(10) 


0.2762(1) 
0.3006(3) 
0.15(3) 


7.88(1)% 
10.96(2) % 
15(1) % 



and (3 states and the bound deuteron state (which we 
denote by a subscript 7 in this section), 

POO 

= (7 2 + fc 2 ) / dr u 7 (r)u fe (r) + w y (r)w k (r) 
Jo 1 J 

= [u'^Uk - Uju' k + w'^Wk - Wjw' k ] ^ (38) 
Using the short distance solution, Eq. I|14|) . we get 

CU,i(k)CU(7)sinfo>(7) - <Pi(k)] =0 , i = a,p (39) 
Which yields 

^(7) = <p a (k) = Mk) (40) 

Thus, the short distance phases <p a (k) and <pp{k) of the 
3 S± — 3 Di channel wave functions in the OPE potential 
at short distances are all determined by deuteron prop- 
erties. This means, in particular that the low energy 
parameters and scattering phase shifts are uniquely de- 
termined by the deuteron binding energy and the OPE 
potential parameters 5 . 

The previous argument can also be implemented if we 
have a short distance cut-off at r = a, the orthogonality 
relation of Eq. (|38|) transforms into the condition 

u' 7 (a)u k: i(a) + w^(a)w fc! i(a) (41) 

= u ii a ) u 'kA a ) + w ~/( a Wk/ a ) 

i = a,0 (42) 

Thus, if we impose the same condition on both solutions, 
Eq. H42|) cannot be satisfied unless they are related at 
the boundary. For instance for the condition BC6 of 
Eqs. (|15fl . we get the two relations 



u 'kA a ) = ^ w 'kA a ) 



a, (3 



(43) 



5 This property does not hold for other triplet channels with higher 
partial waves, because there are no bound states in those chan- 
nels. Nevertheless, it is also true that there is only one indepen- 
dent parameter. This means in practice that one can use one 
scattering length out of the three to predict the phase shifts also 
in other partial waves. 



The orthogonality relation corresponding to boundary 
conditions of the form of Eq. I|43|l implies then the or- 
thogonality constraint 



J k.i 



(a) 



w'Ja) 



Uk,i{a)V2 + Wk,i{a) u 7 (a)V2 + w 7 (a) 



(44) 



Which is the analog finite cut-off condition of Eq. I|40|l . 
The remaining conditions in Eqs. i|15|l generate analogous 
orthogonality constraints. 

The results for the 3 Si — 3 Di channel phase shifts using 
these conditions are presented in Fig. [S] The description 
is rather satisfactory and it seems to work, as one might 
expect, up to the vicinity of the CM momentum which 
magnitude coincides with the two pion exchange left cut 
k = i m. 



B. Low energy parameters 

In the low energy limit one has 



Si 

52 



-a k , 

— Ct2k 5 . 

— k 2 

a 



(45) 



so that the zero energy the wave functions behave asymp- 
totically 



uo,a(r) -> 


1-1- 




wo,a(r) -> 


3a 02 


9 5 




UQ,p(r) -> 


r 




ao ' 


w o,/3(r) = 


3a 2 




ao2r 2 



15ai 



1)2 



(46) 



Using these zero energy solutions one can determine the 
effective range. The 3 5i effective range parameter is 
given by 



ro = 2 



1 - 



r 



- u a (r) 2 - w a {r) 2 



dr . 



(47) 
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FIG. 5: Eigen Phase shifts for the OPE potential as a function of the CM np momentum in the triplet 3 <Si — 3 Di channel 
compared to the Nijmegen results piftl ]. The regular scattering wave functions are orthogonal to the regular deuteron bound 
state wave functions constructed from the OPE with 7 = 0.231605fm _1 , m = 138.03MeV and g^NN = 13.1083. 



In the zero energy case, the vanishing of the diverging 
exponentials at the origin imposes a condition on the a 
and (3 states which generate a correlation between a , 
«02 and a 2 . Using the superposition principle of bound- 
ary conditions we may write the solutions in such a way 
that 

UqJt) = ui{r) —u 2 (r) + ^^-u 3 (r) 

tto.aW = Wr) w 2 (r) + ^^-w 3 (r) 

a a 

up{r) = —u 2 (r) + ^^u 3 (r)-—^ — uAr) 
a a 02 15a 02 

w 0(r) = — w 2 {r) H w 3 (r) - — w 4 (r) 

a a 02 15a 02 

where the functions ^1,2,3,4 and ^1,2,3,4 are independent 
on ao, olq 2 and 012 and fulfill suitable boundary condi- 
tions. As a consequence we get a linear correlation be- 
tween 1/aoj otQ 2 /aQ and also a linear correlation between 
a 2 /ao2 and l/ao2- This means in turn that according to 
the OPE potential both a^ 2 and a 2 depends linearly with 
ao- Numerically we get the following correlations, 

a 02 = 0.963571370240 a - 3.467616391389 
a 2 = 3.467616391389— + 5.080264230656 (48) 

These relations are cut-off independent and unique con- 
sequences of the OPE potential. On the other hand, the 
orthogonality between the bound state and the scattering 
state yields 

a = 1.037805911852 a 2 + 3.598712446758 (a) 
a 02 = 0.288382561043 a a 2 - 1.465059639612 a (fi) 

(49) 

The provided high accuracy is indeed needed. The four 
equations, Eq. I|48|l and Eq. (|49|l . overdetermine the val- 
ues of the three scattering lengths and could be solved 



in triplets yielding four different solutions. Actually, 
there are only two independent solutions which differ- 
ences are compatible within our numerical uncertainties. 
The scattering lengths and effective range are presented 
111 TableHJand compared to their perturbative value ( see 
AppendixE] and to the high quality Nijmegen potential 
models [23] 6 - As we see, the agreement with the high 
quality potentials is at the few percent level. Perturba- 
tion theory does not account for most of the contribution 
to the effective range since the orthogonality constraints 
preclude a short distance contribution to r and also to 
the deuteron matter radius r m . This means in practice 
that the counterterm named C 2 in Rcfs. 16, 36| must 
vanish (See Appendix lAl for a detailed discussion). The 
dependence of the scattering lengths ao, ao2 and a2 on 
the deuteron wave number 7 can be seen in Fig. EI where 
7 dependent generalizations of the correlations, Eq. 148(1 
and Eq. gHl bold. 

The phase shifts look very similar to previous work plf 
using an energy expansion of a square well potential as 
a counter-term and adjusting the depth of the two low- 
est orders to reproduce the 3 5i scattering length ao and 
effective range vq as independent parameters or our vari- 
able phase approach with non-trivial initial conditions in 
Ref. [35| where the full coupled channel S-matrix was tai- 
lored to reproduce the effective range expansion to any 
order treating all parameters as independent. Roughly 
speaking, both approaches could be mapped to an en- 
ergy dependent boundary condition with no a priori or- 
thogonality constraints 7 . The fact that the orthogonality 



6 The values of ao and vq have been determined in Ref. 1371 . 
whereas 002 and «2 have been determined by us in Ref. l35l . See 
also Ref. l3tH for a extensive determination in all partial waves. 

7 The main difference in this regard has to do with the multi- 
valuation problem of the potential counter-term in Ref. I2T| typ- 
ical of inverse scattering problems. The approach of Ref. |3~1 
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FIG. 6: The dependence of the 3 Si — 3 Di scattering lengths ao > a 02 and 02 on the deuteron wave number 7 (in fm 1 ). The 
point represents the experimental values. 



constrained boundary condition generates the bulk of the 
low energy threshold parameters with only one parame- 
ter naturally explains the similarity between the present 
phase shifts and those in previous works [2^ and suggests 
that there is perhaps no need to make the short distance 
boundary condition energy dependent if the short dis- 
tance cut-off is removed. 



VI. SHORT DISTANCE SOLUTIONS AND 
DETERMINATION OF THE COEFFICIENTS 

In this section we determine the coefficients of the OPE 
deuteron wave functions appearing at short distances in 
Eq. (|ll|l . In particular, we compute the energy indepen- 
dent and OPE potential parameters independent short 
distance phase ip. Let us remind that any choice of ip cor- 
responds to a different choice of short distance physics; 
given ip and the OPE potential all deuteron and scat- 
tering properties are uniquely determined. However, the 
leading asymptotic form cannot directly be used to match 
the numerical solution obtained by integrating in the 
large distance solution. On the one hand, if we use cut-off 
approaches to determine the regular solution, there are 
short distance cut-off effects when the distance gets close 
to the cut-off radius. On the other hand, the fact that the 
diverging exponential dominates over the converging one 
provides too weak a signal for the corresponding coeffi- 
cient. To remedy the situation we improve on the short 
distance solution to provide a reliable approximation at 
larger distances (~ lfm) where the diverging exponential 
is less dominant, and look for plateaus in the matching 
radius. It turns out (see below) that one should go at 
eight order in this expansion for a robust determination 
of the short distance coefficients. Actually, one can then 



does not have this problem. 



directly match the short distance improved wave func- 
tions to the numerical solution without no reference to 
cut-offs. We will try the two methods and see that they 
yield to compatible results for the short distance coeffi- 
cients. 

In the limit r — ► the solutions to the coupled equa- 
tions can be written in an expansion of the form 8 

u(r) = «o(^) ai e a °Vff( r ) 
w{r) = w o{^) e»»V7 S ( r ) 

(50) 

with 

71=0 

00 In 

t-^ / r \ n 2 
9ir) = EM*) 

n— 

(51) 

At leading order we get the equations 

MoOg + 16V2wo = 
16v^u + Oo - 16 ) w o = (52) 



This expansion looks similar to a coupled channel WKB expan- 
sion but it is free of some inconveniencies. The applicability 
condition of the coupled channel WKB method would be that 
the de Broglie local wavelength matrix should be a slowly vary- 
ing function of distance, implying in turn three conditions on 
the corresponding local wavelength eigenvalues as well as the 
corresponding WKB mixing angle which need not be necessarily 
satisfied simultaneously, generating conversion mode problems. 
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which have the four non trivial solutions. 



(IA) , a Q 

(2A) , a 

(2R) , a 

(1R) , a„ 



-4i, 

+Ai , 

-4n/2, 
+4V2, 



w 
w 



Up 

up 
V2' 



w = -V2u , 
w = -v2wo ■ 



(53) 

(54) 

(55) 
(56) 



The next to leading order equation becomes compatible 
only if 



the range between 0.6 and 1 fm 9 . Matching directly 
the integrated in solution to the short distance solution 
with a vanishing coefficient of the diverging exponential 
C\b. = 0, we get at the scale 0.7 < r < 0.9fm 



C 2R = -0.47(1), 
Cia = 0.1327(3) 



r, = 0.0263333(1) 
,C 2A = 0.2277(5) 



(59) 



We can instead determine rj — 0.0263332 from the bound- 
ary condition BC6 in Eq. (|15fl at r = 0.2fm and deduce 
the remaining constants yielding 

\Ci R \ < HT 7 , C 2R = -0.47(1) 
C 1A = 0.1327(3), C 2A = 0.2277(5) (60) 



a\ = a 2 — 3/4 



(57) 



For any solution in Eq. I|52|l we may then solve for the 
remaining coefficients. One peculiar feature of this ex- 
pansion is that if one wants to determine the solution 
to a given order, one has to compute the coefficients at 
a higher order. The reason is that strictly speaking a 
truncation of the expansion involves also non-diagonal 
elements, and one has the freedom to choose between 
solving u or w to a given accuracy. The explicit result 
to eight order is presented in Appendix El The general 
short distance solution is written as a linear combination 
of the four independent solutions, 



u(r) 



3/4 



w(r) 



V3 Vi?/ 

C 2R f 2R (r)e- 4 ^^ + V2C 1A fiA(r)e 
V2C 2A f2A(r)e il ^' 

V2C 2R g 2R (r)e- 4 ^^ + C 1A g 1A (r)e 
C 2A92A (r)e 4 ^ 



4i\ 



-4n & 



(58) 



This expansion converges rather fast for each solution up 
to distances of about r ~ 0.6 — 0.9fm. That is about 
what one needs, since that is sufficiently far above the 
cut-off radius a ~ O.lfm. Matching u,w,u' and w' at 
some point in this region we get a linear relation between 
Cir, C 2R , C\Ai C 2A and r\. Actually, we find that the 
signal of the converging exponential is about hundred 
to thousand times that of the diverging exponential in 



The errors have been estimated by varying the matching 
point in the region 0.7fm < r < 0.9fm. Note that al- 
though the coefficient of the diverging exponential C\ R 
is six orders of magnitude larger than the one of the 
converging exponential, the solution through the match- 
ing condition is about eight orders of magnitude smaller 
(may even change sign ). So that the result provides a 
sizeable signal for the converging exponential. With these 
values we show in Fig. [7| the short distance wave func- 
tions compared to the integrated in numerical ones when 
the matching is undertaken at r = 0.8fm. To improve on 
the short distance side we have taken C\ R = 0. The error 
in the region 0.7fm < r < 0.9 never exceeds a 0.01%. We 
have checked that setting the constant C 2R = intro- 
duces a larger deviation from the numerical solution as 
compared to the computed value in the region above lfm. 
Finally, the corresponding short distance angle reads 



tan 



-i C 2A 
C\ A 



-59.7(1) 



(61) 



The discussion in this section explicitly shows that con- 
trary to the findings in Ref. [TlJ the coefficient of the 
converging exponential does not vanish. 



9 For instance at r = 0.8fm we get 

u = 1.9683 -59.4526 t; 

= 1.14054 C XA -52.3866 C lR + 

w = -4.00531 + 159.28 ?7 

= -0.667631 C 1A - 551.52 C 1R 

u' = -6.84126 + 287.992?? 

= -1.34949 Cia +273.671 C 1R - 

w = 16.2726 - 607. 14 rj 

= -0.667631 C 1A - 551.52 C 1R 



1.09959 C 2A - 0.00169556 C 2R , 



1.67549 C 2A + 0.0177228 C 2R 



4.02872 C 2A - 0.00925419 C 2R 



1.67549 O 



2 A 



0.0177228 Co 



where the l.h.s. corresponds to the numerical solution and the 
r.h.s. to the short distance approximation, and the barred co- 
efficients Cia = (Cia + C 2A )/2 and C 2A = (~C 1A + C 2A )/2i 
have been introduced. 
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FIG. 7: The short distance expansion for the deuteron wave 
functions matched to the numerical solution at a distance 
about r = lfm. We take C 1A = 0.1328, C 2 a = 0.2277, C 2 r = 
—0.46 and Cm = 0. The numerical solutions are normalized 
by taking As = 1. 



VII. CONCLUSIONS AND OUTLOOK 

In this paper we have reanalyzed the OPE potential 
in the triplet 3 Si — 3 Di channel both for bound and 
scattering states. Rather than modeling the interaction 
below some finite short distance we have adopted the 
viewpoint of taking the potential seriously down to the 
origin. This must be carefully done and in a way as 
to get rid of any short distance ambiguities. In addi- 
tion, this procedure proves crucial to be able to disentan- 
gle the OPE contribution from other contributions, like 
TPE and higher, electromagnetic effects and relativistic 
corrections to deuteron and NN scattering observables. 
Our analysis is carried out entirely in coordinate space 
where these corrections generate a potentical which is fi- 
nite everywhere except at the origin. Momentum space 
treatments require an additional regularization of the po- 
tential. 

The OPE coupled channel potential is singular at short 
distances and additional conditions need to be specified 
on the wave functions at the origin. Actually, the singu- 
lar eigen-potentials at short distances are attractive and 
repulsive and while in the attractive case a mixed bound- 
ary condition specifies the corresponding short distance 
eigenfunction, in the repulsive case one must impose a 
standard homogeneous boundary condition. This only 
leaves one free parameter, which we have chosen to be 
the deuteron binding energy and which cannot be deter- 
mined from the OPE potential. All remaining deuteron 
observables come out for free. For the scattering states in 
the 3 Si — 3 Di channel, we have demanded orthogonality 
constraints between all states of different energy. This 
condition is actually an additional requirement for sin- 
gular potentials, since the orthogonality relation carries 
information on the peculiar short distance behavior of 



the wave functions, and is not necessarily satisfied. The 
most obvious example where orthogonality constraints 
are violated corresponds to energy dependent potentials 
and energy dependent boundary conditions in coordinate 
space. A less trivial but significant example is the case of 
dimensional regularization in the PDS scheme as a per- 
turbative analysis in coordinate space of both bound and 
scattering states reveals. The power of the orthogonality 
constraints for singular potentials is that all scattering 
properties are then predicted from the OPE potential 
parameters and the deuteron binding energy. 

In our analysis it turns out that the short distance form 
of all wave functions is characterized by some short dis- 
tance constants. We have clarified the role played by the 
exponentially suppressed regular solution by determining 
its non-vanishing value numerically using short distance 
expansions to high order, to explore the region below 
O.lfm, not accessible to standard numerical integration 
methods. Another relevant constant is given by a short 
distance phase <p which plays the role of a fundamental 
dimensionless constant in the OPE problem. It does not 
depend on the energy nor on the OPE parameters, but it 
is related to the form of the OPE potential in the chiral 
limit. The closeness of this phase to 7r/3 is mysterious 
and suggestive and requires further investigation. 

It is remarkable that indeed the bulk of the experi- 
mental results both for the bound state as well the scat- 
tering observables are accounted for at the 2 — 3% level 
by the OPE potential taken from zero to infinity. We 
interpret this success as a confirmation on the validity 
of our choice of regular solutions and the use of orthog- 
onality constraints. The discrepancies can legitimately 
be attributed to other effects such as TPE, electromag- 
netic and relativistic corrections. Many of the methods 
and results obtained in this paper can be generalized in 
a straightforward manner to take these effects into ac- 
count and to the study of higher partial waves without 
any substantial modifications. In particular, the num- 
ber of independent constants in a given channel depends 
on the short distance behavior of the long range poten- 
tial. The bonus of such a program would be the complete 
elimination of short distance ambiguities in the study of 
the NN interaction with known long distance forces as 
determined by chiral symmetry. In our view this an in- 
dispensable prerequisite to asses the relevance of chiral 
symmetry in nuclear physics in a model independent way. 
The systematic study of these effects will be reported 
elsewhere p4| . 
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APPENDIX A: PERTURB ATIVE SOLUTIONS 
1. Bound state 

In this appendix we solve the coupled deuteron equa- 
tions, Eq. (|2J) in standard perturbation theory for the 
fixed negative energy bound state. A somewhat re- 
lated approach looking for the equivalence with the PDS 
scheme of 36] in the one-channel positive energy case can 
be looked up in Ref. jijj. The problem of orthogonality 
was not discussed. The requirement of normalizability 
of the deuteron state requires the D wave component to 
vanish. Thus, at lowest order we have the normalizable 
solutions, 



4 0) 

4 0) 



('■) 



(Al) 



At first order we have to solve the equations 



4 1 )"(r)+ 7 2 M «(r) = -U s (r)e-i r . 



(r) 



(r) = -U sd {r)e 



(A2) 



Using the regular and irregular solutions at the origin 
sinh(7r) 





:(r) = 




(r) = 


xXxrrej 


;(r) = 




(r) = 


get 





= 2- 



2 1 



( 7r )5 



sinh^r) 



6 



cosh('yr) 



-7r 



3 

7r 



(A3) 
(A4) 



,.(!) 



W 



(1) 



(r) 
(r) 



G s (r,/)t4(/)e- 7r V 
G d (r,r')[/ sd (r')e-^'dr' 



(A5) 
(A6) 



where G s and Gd are the corresponding Green functions. 
Explicit calculation yields 



u«(r) 



_ 7r m i?r(0, mr + 2 r 7) — T(0, to r) 2 to 2 R Ei(— m r — 2 r 7) sinh(r 7) 



1 + 



3 7 
3 



3 7 



i 2 i? (3 to 2 - 4 7 2 ) T(0, m r + 2 r 7) - T(0, to r) 
6\/27 3 

i? (-6m 3 7 + 6m 2 7 2 + 47 4 + (3m 4 - 4to 2 7 2 ) log(l + ^)) 

6V2 7 3 

e -m r-2 r 7 ^ + 6 m r + to 2 r 2 - to 3 r 3 + 4 r 7 + 4 to r 2 7 + 2 to 2 r 3 7) 

2V2V7 3 

i? e _m r (6 + 6 to r + to 2 r 2 — m 3 r 3 — 8 r 7 — 8 m r 2 7 + 4 r 2 7 2 + 4 to r 3 7 2 ) 



2 V2 
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1 rsj£y£ J ryj> 1 
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2 V2 
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y4 ^y3 



6V2 



(A7) 



(A8) 



where T(0, z) and Ei(z) are the standard incomplete respectively 
Gamma function and the Exponential integral function 



r(o,z) 

Ei(z) 



dt- 



t 



PI dt- 
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At asymptotically large distances we have 



where 



u«(r) 



Cpcrt 
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Cpcrt^ 



??porte 



U s (r)u Icg (r)e ir dr 



U sd {r)w rcg (r)e lr dr 



Explicit calculation yields 
Rm 2 



-•pert 
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log 1 



27 
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6\/2mj 3 



47 4 + 6m 2 j 2 — 6m 3 7 



27 



(3m 4 - V]') log 1 + — 
\ m 

32V2R 2 2V2i?7 3 



45m 



3m 2 



(A16) 



The numerical value we get is ry per t = 0.0510 almost twice 
the exact OPE result. Taking this perturbative value 
for 77 we show in Fig. [8] the perturbative deuteron wave 
functions as compared to the exact ones. 

Unfortunately, if one wants to improve on this first or- 
der calculation going to second order perturbation theory 
there is a problem since the behavior of the perturbative 
wave functions at short distances is given by 



Rm 2 ( 2 7 
— — log 1 + — 

37 \ m 



' r 3 



(A17) 
(A18) 



making the wave function non normalizable, unlike the 
exact regular wave function. This divergence at short dis- 
tances actually precludes going to higher orders in per- 
turbation theory. 



The normalization at first order is given by 

1 r°° 

-72 = / (e-^ r + 2u^(r)e-^) 



(A19) 



and hence 
A s 
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= 1 - 



= 1 - 



2i?m 2 Rm 2 
3(m + 2 7 ) + ^7" 
2i? 7 16i? 7 2 



log 1 + - 



3 9m 
The deuteron matter radius is given by 



(A20) 



r 2 = - A 2 
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r 2 (e- 2 ^ r + 2u i .p{r)e"> r ) (A21) 
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FIG. 8: Perturbative deuteron wave functions compared to 
the exact ones as a function of the distance (in fm ). The 
exact ones are normalized by taking As = 1. 



and hence to first order one has 

2 1 m 2 R(3m + IO7) 



m,pert g 7 2 t l8 7 (m + 2 7 ) 3 
yielding in the weak binding regime 

/- 2R 1 16i? 7 2 

V87 r m)Pert = 1 



(A22) 



(A23) 



3 9m 

Finally, the quadrupole moment at first order is given by 
V2 



'pert 



10./0 ■ 1 



r 2 w^{r)er^ r dr (A24) 



The integral can be evaluated to give 

8R(4m 2 + 97m + 67 2 ) 



Q 



pert 



45 (m + 27) 3 
32i? 8i?7 128i?7 2 
45m 3m 2 15m 3 



+ . . . (A25) 



yielding Q P ort = 0.4555fm 2 . 

Our perturbative expressions for As, r m and Q coin- 
cide with those of Kaplan Savage and Wise provided 
one takes in their expression for r m the renormalization 
scale in the PDS scheme to be /1 = 7, instead of taking 
fx = m as they do or else taking the Ci counter-term 
identically equal to zero. Actually, the Ci counter-term 
can be mapped into a short distance contribution to the 
effective range parameter ro in the 3 Si channel. The 
value of r\ was not given in that reference but can be de- 
duced from the off-diagonal scattering amplitude in the 
3 S*i — 3 Z?i channel given in their previous work |36| by 
evaluating the residue at the deuteron pole. The result 
also agrees with the calculation presented here. 
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2. Low energy parameters 

To check the identification 6*2 = further let us com- 
pute the S-wave effective range r . For our purposes it 
is sufficient to analyze the zero energy scattering state. 
The lowest order solution is given by an a state 



4%) = 



1 



r 



(A26) 



At zeroth order in the OPE coupling the orthogonality 
constraint yields 







7 



+ 7 

a 



(A27) 



which yields the scattering length to lowest order, 



a, 



(0) 



1 

7 



(A28) 



At first order we use the regular solution w rog (r) = r 
and the irregular solution Ui Iles (r) = (1 — r/oto) and get 
similarly to the bound state case the first order correction 
to the a state, 

(1) / s 2i? _ mr , , 2m 2 rR 

u o,L( r ) = 1; — e (1-aom) Ei(-mr) 

dap o 



(!) 1 \ 



R 



15 y/2 ao m 2 r 2 



120 - 64 a m + 120 m r - 34 a m 2 r 



+ 40 m 2 r 2 — 2 oto m r 2 + ao m 4 r 3 — ao m r 4 
R (-120 + 64 a m + a m 6 r 5 T(0, m r)) 



15 a/2 «o m 2 r 2 



(A29) 



Note that asymptotically the first order correction to 
the S-wave vanishes exponentially and hence cannot con- 
tribute to the scattering length. On the other hand, the 
orthogonality relation to first order reads 







dr 



^uZ + u^l+u^l] (A30) 



and after computing the integrals one gets 



1 R 
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Solving perturbatively for «o we get at first order 



1 2m 2 R 



C*0,pcrt — — 
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3 7 2 



7(7 — 2m) 



lo;. ( 1 + 

m 



Numerically one gets a , P ert = (4.3177 + 0.2912 + . . . )fm 
to be compared with the full OPE result ao — 5.34 and 
the experimental value «o = 5.42fm. The E\ scattering 
length «02 can be read off from the D-wave, using the 
asymptotic condition in Eq. 146|) 



ao2 = 



4^2^(15 - 8a m) 

45m 2 
4\/2i?(157- 8m) 
45 7 m 2 



(A33) 



in the second line we have substituted the perturbative 
relation ao = 1/7 + O(R). Note the linear correlation 
a 02 = 1.54 99ao — 4.1530 to be compared with the ex- 
act OPE relation in Eq. I)48[l. The numerical value one 
gets for the first and second lines taking the experimen- 
tal values of ao — 5.42 and 7 are ao2 = 4.24fm 3 and 
ao2 = 2.53fm 3 respectively to be compared with the ex- 
perimental aw = 1.64fm 3 . In the weak binding limit one 
obtains 



7 expert 



1 



2i? 7 16i?7 2 
9m 



(A34) 



In this limit we have the perturbative linear correlation 
between the scattering length and the deuteron matter 
radius 

a 



2V2 



C( 7 3 ,i? 2 



(A35) 



which yields the value r m = 1.92 for the experimental 
scattering length ao = 5.42fm. The linear correlation 
was established empirically with realistic potentials in 
Ref. Ellil. 

To first order the effective range in the 3 Si eigen- 
channel is given by 



fa 



yielding 



?~0,pert 



dr Uol{r)u ( ]l(r) 



4R(3m 2 - 87m + 67 2 ) 
9m^ 

1.4369 - 5.47897 + 5.87587 2 
0.4831fm 



(A36) 



(A37) 



a result much smaller than the full OPE result (1.64 fm) 
and the experimental number (1.75 fm). Again, our re- 
sult corresponds to a theory where the short distance 
contribution to the effective range vanishes, i.e. C2 = 0. 
A non vanishing value of C2 was needed to fit the experi- 
mental values of both the matter radius and the effective 
range. Our calculation shows that the scheme developed 
in Refs. 0] and Ref. j3|| does not fulfill perturbatively 
the orthogonality constraints. 



APPENDIX B: SHORT DISTANCE EXPANSION 



(A32) For the f(r) function we get (we use x = r/R), 
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APPENDIX C: LOCAL DIAGONALIZATION 
AND PERTURB ATIVE MIXING 



One of the puzzles one encounters in the description 
of the deuteron with the OPE potential is that while the 
dimensionless D/S ratio parameter is rather small at long 
distances w/u — > 77 = 0.0256, it actually comes from a 



strong mixing at short distances where w/u — > l/y/2 ~ 
0.707. Actually, the analog question for scattering states 
is that there seems to be a natural hierarchy for the phase 
shifts in the 3 Si — 3 D\ channel, namely 5^si 3> S3D1 3> £1 
even though the threshold behavior of the D-wave is more 
suppressed than that of the mixing angle. The question 
is whether one can think of an expansion in terms of the 
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■q parameter. There are two obvious situations where the 
mixing does not occur. One is the absence of tensor force. 
In the OPE potential that would also eliminate the D — 
wave. Another situation is dropping the mixing terms in 
the OPE potential, which is questionable since they are 
actually larger than the diagonal terms. It is possible, 
however, to write the equations in a form that the mixing 
is manifestly small at all distances. To this end we make 
a local rotation of the deuteron wave functions 




cos 9(r) sin#(r) 
-sin#(r) cos9(r 




(CI) 



in such a way as to diagonalize the potential we have 
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The deuteron equations for the OPE potentials read 
after the local rotation 



<(r) + [U A (r) + 8'(rf] u A (r) + 7 2 u A (r) = [29' (r)u' R (r) + 9"(r)u R (r)} , 
-u£(r) + [U R (r) + 9'{rf] u R {r) + 7 2 Ufl (r) = - [29' (r)u' A (r) + 9"(r)u A (r)} , 
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In the coupled channel space these equations can be vi- 
sualized as a particle with spin in the presence of a gauge 
potential 9'{r). At long distances we have the expan- 



sions 
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Note that in the locally rotated basis the mixing is related 
to the derivative of the mixing angle, 9' which is small 
at all distances (See Fig. EI) • Actually, at asymptotically 
large distances we have 

u A (r) — > u(r) u R (r) — * w(r) (CIO) 

If we neglect the mixing term in Eq. I|C3|I the equations 
decouple and, actually, there is no non-trivial solution 
for the repulsive eigen-channel, since the energy is fixed 
arbitrarily. Hence in the absence of mixing we have u R — 
0. At this level of approximation we then get 
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In Fig. EH we show the solutions of the decoupled equa- 
tions compared to the exact ones. As we see, the differ- 
ence in the wave functions and hence the D/S mixing is 
indeed small. Note that this is not the same as to neglect 
the tensor force. The results for the deuteron observables 
are presented in Table II I II As we see, the quality of the 
zcroth r\ approximation is rather good. 

Actually, we can check a 'posteriori that the mixing is 
indeed small for the zeroth order solutions. The inhomo- 
geneous term at short distances behaves as 



29'(r)u' A (r) 
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which compared to the remaining terms in Eq. (|C3|I can 
indeed be considered small. Under these circumstances, 
the mixing can then be included perturbatively, yielding 



u R {r) = / dr'G R (r, r>) [29'(r')u' A (r') + 9" (r')u A (r')} 



(C14) 

with Gr[t, r') the Green function of the homogeneous 
equation in the repulsive eigen-channel, 

Git(r,r') = w rcg (r)uwg (r')9(r - r) 

+ w re g(rVirrcg(r)6>(r-r') (C15) 

with Wronskian equal to unity and w reg (r) and Wi Ileg (r) 
the regular solution and irregular solutions at the origin 
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FIG. 9: (Left) Reduced eigen potentials (in fm 2 ) [/^(f) and Un(r) as a function of the distance r (in fm). (Right) Local 
mixing angle 9(r) as a function of the distance r (in fm). 
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FIG. 10: The rotated eigenfunctions at zeroth order compared to the exact ones. They correspond to take r\ = 0. 



respectively. Asymptotically one has, 
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To get in practice the coefficient C we start with C = 1 at 
short distances and build the ratio to the asymptotic form 
at a sufficiently large distance. With these conditions the 
solution un{r) at large distances behaves as 



ur(t) 



■qe 



(r — oo) 

(C17) 



drw leg (r) [26'{r)u' A {r) + 6"(r)u A (r)} (C18) 



At short distances we get, from Eq. I|C14|) and using the 
asymptotic forms of Eq. I|C16(1 and Eq. I|C9() . the result 

/ r \ 7/4 



Ufl(r) - CC A (^-J cos [Ay/R/r + a) (C19) 

in agreement with the leading short distance behavior of 
the full solution for the combination u(r) — \f2w(r) (see 
Sect. IVIII . The perturbative value for the asymptotic 
D/S ratio we get is 

r, pm . t = 0.0261 (C20) 
quite close to the OPE exact one, 77ope = 0.0263. 



APPENDIX D: LONG DISTANCE SOLUTIONS 

As a complement to the perturbative treatment of Ap- 
pendix ^ we analyze the bound solutions at long dis- 



21 



TABLE III: Deuteron properties for the OPE potential. We use the non-relativistic relation 7 = ^2fi np B with B — 2.224575(9). 
We compare the ^-expansion at leading order (LO), with standard perturbation theory at Next to leading order (NLO) and 
the exact OPE result. 





7(fm- 1 


V 


A s (fm- 1 / 2 ) 


r m (fm) 


Qd(fm 2 ) 


Pd 


OPE-?? (LO) 


Input 





0.8752 


1.9423 


0.1321 


6% 


OPE-pert (NLO) 


Input 


0.051 


0.7373 


1.6429 


0.4555 





OPE-exact 


Input 


0.02633 


0.8681(1) 


1.9351(5) 


0.2762(1) 


7.88(1)% 


Nijmll 


Input 


0.0253(2) 


0.8845(8) 


1.968(1) 


0.271(1) 


5.67(4)% 


Exp. (non-rel.) 


0.231605 


0.0256(4) 


0.8846(9) 


1.971(6) 


0.2859(3) 


5.67(4)% 



tances. The asymptotic deuteron wave functions for the The first order solution can be evaluated analytically, 
OPE potential can be written in the form yielding 



ulr) = e 



-77" 



5> 



r e 



-kmr 



w(r) = r\e ir H h 



3 3 



7r (7^)' 



L k 



(Dl) 



(D2) 



Fi(r) 



Gi(r) 



Re~ mr (m 2 r 2 - 2 (1 + rj) - 2mr (1 + 7-7)) rj 
V2 r 3 7 2 

m 2 R (3v / 2m 2 ?7 + 7 2 (4 — 4 V2 77)) Ei(-mr) 
12^ 

e 2ri m 2 R (-3\/2?7i 2 77 + 47 2 (-1 + V2rf)) Ei(-mr-2r 7 ) 

12^- 

m 2 R (-47 2 (\/2 — 2 77) +3m 2 (V2-rj)) Ei(-mr) 
12 7 3 7) 

Re- mr (3m 3 r - rj) +m 2 r 1 (-3 + r 7 ) (v^-?/) + 2 7 2 (l + mr) (2j7 + r 7 (-v^ + t?))) 

2 r7 2 (3 + 3r7 + r 2 7 2 ) rj 
e 2r ~*m 2 R (3-3r7 + r 2 7 2 ) (4 7 2 ( v / 2-2t ? ) -3m 2 (v^-r?)) Ei(-mr-2r 7 ) 

12 7 3 (3 + 3 r 7 + r 2 7 2 ) 77 



(D3) 



(D4) 



The second order can also be evaluated but the expres- 
sion is too long to be presented here. In Fig. llll we present 
the first order solutions compared to the exact ones. The 
perturbative solutions of Appendix [S] are obtained from 
the requirement that the S-wave u, becomes normalizable 
when extended down to the origin. This can only hap- 



pen in the D / S asymptotic ratio, 77 takes the value given 
by Eq. (|A16(I . This illustrates the fact that perturbation 
theory can always be applied at long distances but fails 
at short distances. 

Note that here one treats the coupling constant R and 
the mixing parameter 77 as independent variables. 
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